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Abstract 

> 

■ The infinite-dimensional mechanics of fluids and plasmas can be formulated as "non- 

canonical" Hamiltonian systems on a phase space of Eulerian variables. Singularities of 
| the Poisson bracket operator produce singular Casimir elements that foliate the phase 

^ ■ space, imposing topological constraints on the dynamics. Here we proffer a physical 

t^> interpretation of Casimir elements as adiabatic invariants — upon coarse graining mi- 

croscopic angle variables, we obtain a macroscopic hierarchy on which the separated 
action variables become adiabatic invariants. On reflection, a Casimir element may be 
unfrozen by recovering a corresponding angle variable; such an increase in the number 
^ \ of degrees of freedom is, then, formulated as a singular perturbation. As an example, we 

propose a canonization of the resonant-singularity of the Poisson bracket operator of 
the linearized magnetohydrodynamics equations, by which the ideal obstacle (resonant 
Casimir element) constraining the dynamics is unfrozen, giving rise to a tearing-mode 
instability. 



1 Introduction 

Whereas canonical Hamiltonian mechanics is described by a Poisson bracket operator (field 

tensor) that has a full rank on a symplectic manifold, general noncanonical Hamiltonian 

mechanics is endowed with a Poisson bracket operator (henceforth Poisson operator) that 

may have a nontrivial kernel; the corresponding Poisson manifold may then be split into 
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some local symplectic leaves (Lie-Darboux theorem). A Casimir element foliates the Poisson 
manifold (with the gradient of a Casimir element belonging to the kernel of the Poisson 
operator) Consequently, an orbit is constrained to a leaf (level set) of a Casimir element, 
i.e., a Casimir element is a constant of motion. The constancy of a Casimir element is inde- 
pendent of the Hamiltonian (whereas a usual constant of motion pertains to some symmetry 
of a Hamiltonian), and it is due to a singularity of the Poisson operator. Here we proffer 
an interpretation: "a Casimir element is an adiabatic invariant that is separated from a 
microscopic angle variable by coarse graining" — a Casimir leaf is then a macroscopic hier- 
archy. On reflection, a Casimir invariant may be unfrozen by recovering a corresponding 
angle variable. Such an increase in the number of degrees of freedom is, then, formulated as 
a singular perturbation (cf. [25J). 

We will cast the theory of tearing-mode instabilities into a perspective of an unfrozen 
Casimir leaf. The dynamics of an ideal plasma obeys topological constraints on the co-moving 
magnetic fieldjj which prevents the onset of various potential instabilities. A "tearing mode" 
is a mode of deformation that brings about a change in the topology of magnetic surfaces 
(the integral surfaces of magnetic field lines), which is, therefore, forbidden to occur in the 
ideal dynamics [7J [22]. In [27], a tearing mode is formulated as an equilibrium point on a 
helical-flux-Casimir leaf. As long as the helical-flux Casimir is constrained, the tearing-mode 
cannot grow. We will formulate a "perturbation" that allows the system to change the helical 
flux, as well as absorb (dissipate) the energy. An unstable tearing mode has negative-energy, 
with respect to a fiducial "helical" bifurcated equilibrium; hence, a tearing-mode instability 
is a process of moving across Casimir leaves toward a lower-energy equilibrium point with a 
different (helical) topology. 

In the next section, we begin by reviewing some aspects of the basic framework of Hamil- 



1 A Casimir element C is a member of the center of the Poisson algebra, i.e., [C,G] = for all G. For 
finite-dimensional systems, one may regard phase space as the dual of the Lie algebra and Casimir leaves 
as coadjoint orbits. Unfortunately, in infinite-dimensions, there arc functional analysis challenges that limit 
this interpretation (see, e.g., [10]). Note, a Casimir leaf is not necessarily a symplectic leaf, because the 
kernel of a Poisson operator may not be fully integrable when the Poisson operator has singularities (see, 

e.g., mm- 

2 The core of the topological constraint is epitomized by Kelvin's circulation law: in a barotropic fluid, the 
Lie derivative of fluid momentum is an exact differential 1-form; thus, the circulation (co-moving loop-integral 
of the fluid momentum, or the surface- integral of the vorticity) is conserved. In a plasma (charged fluid), 
the momentum combines with the electromagnetic potential (gauge field), and the corresponding canonical 
vorticity is the combination of the fluid vorticity and magnetic field. In the magnetohydrodynamic model, 
the canonical vorticity of the electron fluid is approximated by the magnetic field, neglecting the electron 
mass. The circulation law for the magnetic field implies that the magnetic flux on each fluid element is 
conserved, the so-called Alfven's law. 
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tonian mechanics. In Sec.[3j we then consider an example of magnetized particles, in order to 
establish a connection between adiabatic invariants and Casimir elements. In SeclU we will 
formulate a systematic method of canonization by adding "angle variables" that result in 
the unfreezing of Casimir elements. As is now well known, the infinite-dimensional mechan- 
ics of a plasma can be formulated as a noncanonical Hamiltonian system on a phase space 
of Eulerian variables (see e.g. [IS]). After a short review of the Hamiltonian formalism of 
magnetohydrodynamics (MHD) and its application to the tearing-mode theory (Sec. l5.lj) . we 
will formulate a (formal) singular perturbation that gives rise to a tearing-mode instability, 
and finally discuss its physical implications (Sec. l5.2|) . 

2 Preliminaries: noncanonical Hamiltonian systems and 
Casimir invariants 

We denote by z = (q 1 ,--- ,g m ,p 1 ,--- ,p m ) the state vector, a point in an affine space 
X = M 2m (to be called phase space) u A canonical Hamiltonian system is endowed with 
a Hamiltonian H(z) (a real function on the phase space X) and a 2m x 2m antisymmetric 
regular matrix 



where I m and m are the m-dimensional identity and nullity, respectively. (In what follows, 
we will write just / or without specifying the dimension, especially when we consider 
an infinite-dimensional space.) We call J c the canonical Poisson operator (matrix). The 
equations of motion (Hamilton's equations) are written as 



whence an equilibrium point is seen to satisfy d z H(z) = 0. Defining a Poisson bracket by 



the rate of change of an observable f(z) is determined by 

!/ = [/,*]. 

One may generalize the Poisson operator J to be a function J(z) of an arbitrary dimension 

n x n (here we assume a finite n, while we will consider infinite-dimensional systems later). 

3 Usually, phase space is identified as a cotangent bundle T*M of a smooth manifold AI of dimension m, 
on which a symplectic 2-form oj = (1/2) J c ^idz k Adz £ (the vorticity of a canonical 1-form) defines symplectic 
geometry. 






[a,b] := (d Zi a)Jij(d Z:j b), 
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A noncanonical Hamiltonian system allows J(z) to be singular, i.e., Rank J(z) may be less 
than n and can change as a function of z (while the corresponding Poisson bracket must 
satisfy Jacobi's identity). The equations of motion are, then, 

±z = J(z)d z H(z). (2) 

A Casimir element C(z) is a solution to a partial differential equation (PDE) 

J(z)d z C(z) = 0, (3) 

which implies that [C, F] = for every F. Therefore, C is a constant of motion (dC/dt = 
[C, H] — for any Hamiltonian H). 

Obviously, if Rank J(z) = n (the dimension of the phase space), ([3]) has only the trivial 
solution (C = constant). If the dimension v of Ker(J(z)) does not change, the solution 
of d2J) may be constructed by "integrating" the elements of Ker(J(z)) — then the Casimir 
leaves are symplectic manifolds. This expectation turns out to be true provided the Poisson 
bracket satisfies the Jacobi identity and m — v is an even number (Lie-Darboux theorem). 
However, the point where the rank of J(z) changes is a singularity of the PDE ([3]) [16], from 
which singular Casimir elements are generated [26] . 

When we have a Casimir element C(z) in a noncanonical Hamiltonian system, a trans- 
formation of the Hamiltonian H<z) such as 

H(Z) ^ H^z) = H(z) - (iC(z) (4) 

(with an arbitrary real constant fi) does not change the dynamics. In fact, Hamilton's 
equations (J2]) are invariant under this transformation. We call the transformed Hamiltonian 
H^z) an energy- Casimir function [HI tUl HSl E] • 

Interpreting the parameter fjj cLS cL Lagrange multiplier of the equilibrium variational 
principle, H^(z) is the effective Hamiltonian with the constraint that restricts the Casimir 
element C(z) to be a given value (since C(z) is a constant of motion, its value is fixed by its 
initial value). As we will see in some examples, Hamiltonians are rather simple, often being 
"norms" on the phase space. However, an energy-Casimir functional may have a nontrivial 
structure. Geometrically, H^z) is the distribution of H(z) on a Casimir leaf (hyper-surface 
of C(z) = constant). If Casimir leaves are distorted with respect to the energy norm, the 
effective Hamiltonian may have a complex distribution on the leaf, which is, in fact, the 
origin of various interesting structures in noncanonical Hamiltonian systems. 
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3 Foliation by adiabatic invariants 



Here we study an example of noncanonical Hamiltonian mechanics (and creation of inter- 
esting structures on Casimir leaves) in which Casimir elements originate from adiabatic 
invariants. 

The Hamiltonian of a charged particle is the sum of the kinetic energy and the poten- 
tial energy: H = mv 2 /2 + q<f>, where v := (P — qA)/m is the velocity, P is the canonical 
momentum, (0, A) is the electromagnetic 4-potential, m is the particle mass, and q is the 
charge. Needless to say, a magnetic field does not change the value of energy, and the stan- 
dard Boltzmann distribution function is independent to the magnetic field. However, in the 
vicinity of a dipole magnetic field rooted in a star or planet, for example, we often find a 
plasma clump with a rather steep density gradient. In such a situation, so-called inward 
diffusion drives charged particles toward the inner higher- density region, which is seemingly 
opposite to the natural direction of diffusion (normally, diffusion is a process of flattening 
distributions of physical quantities). Creation of such a macroscopic structure can be ex- 
plained only by delineating a fundamental difference between a macroscopic hierarchy and 
basic microscopic mechanics. Since the magnetic field does not cause any change in the 
energy of particles, there is no way to revise the energy in the calculation of the equilibrium 
state. Instead, the problem is solved by finding an appropriate "phase space" (or an ensem- 
ble) on which the Boltzmann distribution is achieved; the identification of an appropriate 
macroscopic phase-space is nothing but the formulation of what we call a "scale hierarchy" . 

Magnetized particles live in an effective phase space that is foliated by adiabatic invariants 
associated with periodic motions of particles. Denoting by V\\ and v± the parallel and 
perpendicular (with respect to the local magnetic field) components of the velocity, we may 
write 

H = -v\ + -v\ + q<j>. (5) 

The velocities are related to the mechanical momentum via p := mv, p» := mvn, and 
p± := mv±. In a strong magnetic field, v± can be decomposed into a small-scale cyclotron 
motion v c and a macroscopic guiding-center drift motion Vd- The periodic cyclotron motion 
v c can be "quantized" to write mv 2 /2 = fiu c (x) in terms of the magnetic moment \i and 
the cyclotron frequency uj c (x); the adiabatic invariant fi and the gyration phase d c := u c t 
constitute an action-angle pair. The macroscopic part of the perpendicular kinetic energy 
is expressed as = (Pg — qip) 2 / (2mr 2 ) , where Pg is the angular momentum in the 6 

direction and r is the radius from the geometric axis. In terms of the canonical-variable 
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set z = (i9 c , fi, C)P||) 0, Pe)-, the Hamiltonian of the guiding center (or, the quasi-particle) 
becomes 

Note that the energy of the cyclotron motion has been quantized in term of the frequency 
u c (x) and the action /i; the gyro-phase i9 c has been coarse grained (integrated to yield 2tt). 

Now, we formulate the "macroscopic hierarchy" on which charged particles create a ther- 
mal equilibrium. The adiabatic invariance of the magnetic moment /j, imposes a topological 
constraint on the motion of particles; it is this constraint that is the root-cause of a macro- 
scopic hierarchy and of structure formation. The Poisson operator on the total (microscopic) 
phase space, spanned by the canonical variables z = ($ c , /x, (,p\\, 0, P s ), is a canonical sym- 
plectic matrix: 

"-(Hi)- *=-(-" (?) 

The equations of motion for the Hamiltonian H c are written as dz J /dt = [z\H c ]. Notice 
that the quantization of the cyclotron motion in H c suppresses change in /i. 

To extract the macroscopic hierarchy, we "separate out" the microscopic variables (i!? c , /x) 
by modifying the Poisson operator as follows: 



Jnc 



The Poisson bracket [F, G] nc := (d z F, J nc d z G) determines the kinematics on the macroscopic 
hierarchy; The corresponding kinetic equation d t f + [H c , f] nc = reproduces the familiar 
drift-kinetic equation. The kernel of J nc makes the Poisson bracket [, ] nc noncanonical [To] . 
Evidently, \i is a Casimir element (more generally C = g(n) with g being any smooth 
function). The level-set of /i, a leaf of the Casimir foliation, identifies what we may call the 
macroscopic hierarchy. 

By applying Liouville's theorem to the Poisson bracket [ , ] nc , the invariant measure on the 
macroscopic hierarchy is d A z = d 6 z/(27rd/i), the the total phase-space measure modulo the 
microscopic measure. The most probable state (statistical equilibrium) on the macroscopic 















Jc 











Jc 



z, 



ensemble maximizes the entropy S = — J f log / d 6 z for a given particle number N — J fd' 
a quasi-particle number M = j fifd 6 z, and an energy E = j H c fd 6 z. Then, the distribution 
function is 

/ = /„:= Z- X e-^ H ^\ (9) 
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where a, /3, and \ogZ — 1 are, respectively the Lagrange multipliers on M, E, and N. In 
this grand-canonical distribution function, a/ (3 is the chemical potential associated with the 
quasi-particlesO The factor e~ afl in f a yields a direct u c dependence of the coordinate-space 
density: 



which demonstrates the creation of a density clump near the dipole magnetic field [28] . 



4 Canonization atop Casimir leaves 

The aim of this section is to formulate a systematic method of "canonization" of a noncanon- 
ical system by embedding the system into a higher- dimension phase space; Casimir elements 
become "adiabatic invariants" associated with a symmetry (at a macro-scale hierarchy) of a 
Hamiltonian. 

4.1 Extension of the phase space and canonization 

Let J be a Poisson matrix on an n- dimensional phase space X = IR n parameterized by 
z = (z%, • • • , z n ). We assume that Ker( J) has a dimension v and n — v is an even number. 
We also assume that Ker( J) is spanned by Casimir invariants C%, ■ ■ ■ , C u , i.e. 



Our mission is to find the "minimum" extension of the phase space and a canonical Poisson 
matrix J by which the Casimir invariants are re-interpreted as adiabatic invariants — an 
appropriate perturbation of the Hamiltonian will then give a near-integrable system in the 
vicinity of the original Casimir leaves. The target phase space must be of dimension n := n+v 
(even number) consisting of Z\, ■ • ■ ,z n and additional • • ■ , "& v . 

Before formulating such a minimum system, we note that we may formally produce a 
"larger" system; the simplest method of extension and canonization is to double the phase 
space: let z x2 := (zi, • • • ,Z n ,Xi,--- ,Xn), and 



We can also derive (|9]) by an energy- Casimir function. With a Casimir element fx, we can transform 
the Hamiltonian as H c i— > H a := H c + a/i (a is an arbitrary constant) without changing the macroscopic 
dynamic. The Boltzmann distribution with respect to H a is equivalent to (j9]). This equivalency was discussed 
in greater generality in [15] . 




(10) 




Ker( J) = {Vd, • • • ,Va}- 





(12) 
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where L(x) is a certain regular n x n matrix. To satisfy the Jacobi identity, L(x) must 
satisfy the Maurer-Cartan equation; see Eq. (292) of [T6] . 



4.2 "Minimum" canonization invoking Casimir invariants 

It is generally difficult to reduce J X 2 to J of dimension h x n; to separate In — h variables 
from z X 2, these variables and the remaining n variables must be independent, implying a 
"separation of variables." 

Our strategy is to use the Casimir foliation (TITT) of the phase space. We first canonize J 
on X/Ker(J). Let 

z = (Ci) ■ • • > Cn-v, Ci, • • • , CV) e M n , 

by which J is transformed into a Darboux standard form: 



J' 



\ 



We can extend J' to an n x n canonical matrix such that 



/ 



V 



Jr. 



The corresponding variables are denoted by 



(13) 



(14) 



Zex — (Cl; ' ' " i Cn-w, C\i - - • j C u , ■ ■ • , f^) G M n . 

An interesting property of this extended, canonized Poisson matrix J ex is that the el- 
ements are independent of the additional variables i9, which is in marked contrast to the 
simple extension J X 2 defined in ( fl2l) . 



5 Application to tearing-mode theory 

In this section, we put the method of unfreezing Casimir elements to the test by studying the 
tearing-mode instability from the perspective of the noncanonical Hamiltonian formalism. 
The system is of infinite dimension, hence the formulation needs an appropriate functional 
analytical setting. Here we invoke a simple incompressible ideal MHD model. 



S 



5.1 Helicity and Beltrami equilibria 

5.1.1 Magnetohydrodynamics (MHD) system 

Let V and B denote the fluid velocity and magnetic field of a plasma. Here we consider 
an incompressible flow, V ■ V = 0, hence both V and B are solenoidal vector fields. The 
governing equations are (in the so-called Alfven units) 

d t V - V x (V x V) = -Vp+ (V x B) x B, 

d t B = V x (V x B). ' ^ ' 

where p denotes the fluid pressure. We consider a three-dimensional bounded domain fl 
surrounded by a perfectly conducting boundary dfl; the boundary conditions are (denoting 
by n the normal trace onto dQ) 

n-V = 0, n -B = 0, (on 90). (16) 

The state vector u = t (V, B) belongs to the phase space X = L%(VL) x L^(f2), where 

L*(fi) := {u e L 2 (n); V • u = 0, n ■ u = 0}, (17) 

which is a closed subspace of L 2 (Q) (we endow the Hilbert space X with the standard 
L 2 inner product (u,v) and the norm ||w||). We denote by V a the projector onto L 2 (Q). 
Defining a Hamiltonian and a Poisson operator by 

H{u) := \{\\V\\ 2 +\\B\\ 2 ). (18) 
, , ._ ( -V a (V*V)x P„(Vxo)xB\ 

JW ^ vx[oxb] o y [ J) 

the MHD system f|T5|) is cast into the following Hamiltonian form: 

d t u = J(u)d u H(u), (20) 

(cf. [T3l E] ) where d u is the gradient (of Lipschitz continuous functionals [1]) in the Hilbert 
space X. Here, we define J{u) on a subdomain of C°°-functions in the phase space X, which 
suffices to find regular equilibrium points (cf. [26] for more precise definitions). 

5.1.2 Beltrami eigenfunctions 

The Poisson operator Jiu) has two independent Casimir elements (denoting by A the vector 
potential of B) 

Ci(u) :— - f A B d 3 x, C 2 {u) := [ V ■ B d 3 x, (21) 
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which, respectively, represent the magnetic helicity and the cross helicity. They impose 
topological constraints on the field lines [12] . The "Beltrami equilibrium" is an equilibrium 
point of the energy- Casimir functional H{u) — / u, 1 Ci(n) — /^C^**)- Here we consider a 
subclass of equilibrium points assuming //2 = 0. Then, V = (invoking \xi ^ 0, we obtain a 
larger set of equilibria with a finite V). The determining equation for B is (denoting \i\ = fi) 

V x B - fiB = 0, (22) 

which reads as an eigenvalue problem of the curl operator [23]. The solution (to be denoted 
by B^) is often called a Taylor relaxed state [2"Ul |2"T] . 

While the Beltrami equation ( 1221) together with the homogeneous boundary conditions 
(FTO]) are seemingly homogeneous equations, there is a "hidden inhomogeneity" when Q is 
multiply connected [then, the boundary conditions ffTol) are insufficient to determine a unique 
solution]. To delineate the "topological inhomogeneity" of the Beltrami equation, we first 
make Q into a simply connected domain Q$ by inserting cuts across each handle of Q: 
Qs '■= ^ \ (U^ =1 £f) (where v is the genus of The fluxes of .B are given by (denoting 
by dcr is the surface element on S^) $^(-B) := J* s .B • dcr, which are constants of motion. 
To separate these fixed degrees of freedom, we invoke the Hodge-Kodaira decomposition 
Ll(n) = L|(fi) © L|(fi), where 

L|(fi) := {u G L 2 (fi); V ■ i* = 0, n ■ u = 0, $ e (u) = (W)}. (23a) 
:= {u G L 2 (fi); V x u = 0, V ■ u = 0, n ■ it = 0}. (23b) 

The dimension of L^(Q), the space of harmonic fields (or cohomologies) , is equal to the 
genus v of f2. We decompose the total B G £ 2 (fi) into the fixed harmonic "vacuum" field 
-E?h -^h(^) (which carries the given fluxes $i, • • • , <& u ) and a residual component B-% driven 
by currents within the plasma volume Q, 

B = B^ + B u , [S s := VxB G B H G (24) 

where denotes the orthogonal projector from L 2 (f2) onto L|(f2). 

Now the Beltrami equation (1221) reads as an inhomogeneous equation (denoting Vx by 
curl) : 

(curl - /i)B s = //B Hj (25) 

where the harmonic field B^ is uniquely determined by the fluxes $i, • • • ,<Ev When B^ 
and are given, we solve (j2"51) for £? s to obtain the Beltrami magnetic field B^ = B% + B^. 
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If -B H = 0, ( 125]) has solutions only for discrete eigenvalues /x £ {Ai, A 2 , • • • } =: cr p (S) of the 
self-adjoint curl operator S defined on the operator domain [23] 



D(S) = H^(n) := {u £ L|(fi) n i?" 1 ^); V x w £ z|(ft)}. 



(26) 



If 7^ 0, (125]) has a nontrivial solution for every fi ^ o~ p (S) [23]. Moreover, if the vector 
potential A H of B H and the eigenfunction u>j of S belonging to an eigenvalue Xj are orthog- 
onal (i.e. (Ajj,ujj) = 0, the inhomogeneous equation (1251) has a solution Gj at ii = Aj even 
with Bh ^ 0. Then, \x = Xj is a bifurcation point of two branches of Beltrami fields, 
with ii > Aj and -Ba^q = Gj + «^ (a G 1), and the latter has a smaller energy for a given 
helicity C\ and B H [27] . 

5.1.3 Linearization near the Beltrami equilibrium and tearing mode 

In the neighborhood of a Beltrami equilibrium, we find an infinite number of Casimir elements 
stemming from the resonant singularity of the Poisson operator, which foliate the phase space 
and separate the bifurcated Beltrami equilibria on a common helicity leaf. 

We linearize the MHD equations. Since the Beltrami equilibrium = t (0,B^ l ) is a 
stationary point of the energy-Casimir functional = H — fiC±, the linearization of Hamil- 
ton's equation is rather simple: Denoting by u = t (V,B) the perturbed state vector, we 
define linearized Hamiltonian and Poisson operators by 



Evidently, is a self-adjoint operator for every fi £ R. The linearized Hamiltonian equation 
reads 



In what follows, we assume /i > 0. Then, the positive side of the spectrum o~ p (S) plays an 
essential role; for fi < 0, we switch to the negative side of cr p (S). Evidently, \i > Ai destroys 
the coercivity of ('H fl u,u) with respect to the norm ||£t|| 2 , violating the sufficient condition 
of stability [21] (see also [91 [H HU UHl EH])- In fact, a perturbation B oc ui\ (the eigenfunction 
corresponding to Ai) yields {T-L^u, u) < 0. However, the negative energy of a perturbation 
B oc wi does not necessarily cause an ideal-MHD instability, since motion including oji may 
be "inhibited" in the Hamiltonian mechanics. 




(27) 



(28) 



d t u = JuHuU. 



(29) 
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Let us see how Casimir elements foliate the phase space of perturbations: Ker(j7 M ) consists 
of two classes of elements: t (v, 0) and *(0, b) with v and b satisfying, respectively, 

V x (Bp xt))=fl, V ■ v = 0, (30a) 
Bp x (V x b) = 0. (30b) 

The Casimir elements are, in terms of such v and b, 

C v (u) := [ V -vd 3 x, C b (u) := [ B -bd 3 x. (31) 



Obviously, we can choose v = Bp and b = Bp. However, far richer solutions stem from the 
singularity of J p. 

Here we concentrate on the "magnetic part" (130bl) . but a similar singular solution v 
can be constructed for the "flow part" fl30al) . The determining equation f!30bl) of b can be 
rewritten as 

V x b = r]Bp (32) 

with some scalar function rj. We have already found a solution b = Bp and 77 = Here we 
seek solutions with non-constant 77. However, 77 is not a free function; the divergence of both 
sides of (|32|) yields 

Bp ■ Vrj = 0, (33) 

which implies that rj is constant along the magnetic field lines. For the integrability of rj, the 
equilibrium field Bp must have integrable field lines; a continuous spatial symmetry guaran- 
tees this. Here we consider a slab geometry, in which we may write Bp = 1 (0, B y (x), B z (x)). 
Denoting b = * (0, b y (x), b z (x)), f!30bj) reads as 

B y d x b y + B z d x b z = 0, (34) 

which may be solved for b y (x), given an arbitrary b z (x). Furthermore, we have singular 
(hyper-function) solutions; let us consider 

b = 1 (0, b y (x), b z (x)) e< k vv+ k * z ). (35) 

Putting b y (x) = ik y $(x) and b z (x) = ik z $(x), (1511) reduces into 

[B y (x)ky + B z (x)k z ]d x §(x) = 0, (36) 

which yields 

= c + c 1 Y(x-x' ! ), (37) 
12 



where cq, C\ are complex constants, and k y , k z and x^ (real constants) are chosen to satisfy 
the resonance condition 

B y (x ] )k y + B z (x^)k z = 0. (38) 

Then, t] = \{k y j B z )e l( - kyV+kzZ ^5(x — x^). ^From (132]) we see that this Dirac 5-function solution 
implies a current sheet on the resonant surface r' : x = x\ Physically, 1^ represents a thin 
layer of ideal-MHD plasma that supports a sheet current, which suppresses the change of 
field-line topology [2]. 

In what follows, we normalize the kernel element b so that ||6|| 2 = (b,b) = 1. The 
singular (hyper- function) solution b of ( 135|) created by the resonance singularity ( )38|) . imposes 
an essential restriction on the range of dynamics; any magnetic perturbation B such that 
(B, b) 7^ is forbidden to change, because 

C b (u) = C b (B):=(B,b) (39) 

is an invariant. We call Cb{B) a "helical-flux Casimir invariant." The equilibrium point of 
the energy-Casimir functional 

F^(u) := l -{U,u,u) - 0C b (u) (40) 

gives the tearing mode. 

Because of the linearity of the determining equation fl35l) . the totality of *(0, b) G Ker(j7/J 
is a linear subspace of the total phase space and it is "integrable" - thus foliates the phase 
space in terms of the Casimir invariants C b {u) = (B, b). In the next subsection, we choose 
the "dominant helical-flux Casimir" that has the common Fourier coefficients with the helical 
mode oji of the bifurcated helical Beltrami equilibrium, and define the "minimum extension" 
that canonizes the corresponding kernel of J^. 

5.2 Tearing-mode instability 
5.2.1 Canonization 

Let *(0, b) G Ker(j7^). We separate a one-dimensional subspace {pb; p G 1R} from the phase 
space L"^(Q) of magnetic perturbations B, and denote by p\\ the orthogonal projection onto 
the remaining space: 

Pl] B:=B-(B,b)b. 

We also denote 

p ± B := (B,b)b = C b (B)b, 
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and decompose B = p\\B + p±B. Writing the state vector as it' = t (V, fp\\B, p±B), and 
denoting /C M := (1 — fiS' 1 ), the Hamiltonian and Poisson operators read 



J' 

f 1 



1 






p\\K>v. 




(41) 









p||Curl(o x B^ 



(curlp|| o) x B h 







(42) 







Notice that the kernel t (0,b) has been separated from the upper left block of the Poisson 
operator. 

Now, we introduce an adjoint variable q to extend the phase space: 



u, 



'y,P\\B,p\\B,q) 



and define 



/ 



fi,ex 





I curl (o x B p 



(curlpy o) x B^ 




V 



(43) 



-1 

1 / 

which is "canonized" by extending the variable p. Since the original system does not include 
q as an variable, we may write 

/ i o ; \ 

piJ^p, 



n 



)i,ex 



\ 



P±Xn o 
J 



(44) 



Evidently, p±B = (B, b)b is invariant, which is originally a Casimir element, but is now an 
invariant due to the symmetry of H^ex with respect to the new variable q. 

5.2.2 Singular perturbation 

Perturbing the Hamiltonian with respect to q, we can break the invariance of p := (B, b) = 
Cb(B). We consider a Hamiltonian 



( 1 





H,EX 





\\K., 



pj-K-f, 



\ 

.Pwfiv 

PlXpi o 

D J 



(45) 



where D is a parameter that is introduced to couple the original system to the external 
variable q. Note that the original energy ('H fl u, u)/2 is no longer an invariant; instead, the 
new total energy ( / H IJ , ) ExU ex ,u ex )/2 is conserved. 
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The induced change in the Casimir element (helical flux, which is now denoted by p) is 
estimated by the canonized block of Hamilton's equations: 

J t P=~ D ^ ^q = (K lt B,b)- (46) 

For B = puj\ (the eigenfunction determining the bifurcated fiducial-energy equilibrium), we 
may estimate 

(K^B.b) = {(l- t jS- 1 )u li b)p=(l-n/X 1 ){ui 1 ,b)p. 

Absorbing the sign of (u>i, b) by p, we assume 7 := (u>i, b) > 0. For simplicity, let us assume 
that D is a constant number. The system (|46p has the Hamiltonian 

This sub-system Hamiltonian describes the coupling of the original (un-perturbed) Hamil- 
tonian system with an "external energy" Dq 2 /2. If this external energy is positive (i.e. 
D > 0), the "internal energy" of the original system may "dissipate" through the coupling. 
The factor (1 — /u/Ai) of the "kinetic energy" part of the Hamiltonian H p may be interpreted 
as an effective (reciprocal) mass of the tearing mode — beyond the bifurcation point \x = Ai, 
the effective mass becomes negative, and the "negative-energy mode" can grow by absorbing 
energy from the positive energy source Dq 2 /2. 

6 Conclusion 

We have described several facets of noncanonical Hamiltonian systems. Namely, that the 
Poisson operator (field tensor) of a noncanonical Hamiltonian system has a nontrivial kernel 
(and thus, a cokernel) that foliates the phase space (Poisson manifold), imposing topological 
constraints on the dynamics. The Hamiltonian (energy) of a weakly-coupled macroscopic 
system (such as a normal fluid or a plasma) is usually rather simple — it being a convex 
functional (typically a quadratic form) by which one can define an energy norm on the 
phase space. However, an "effective energy" may have a considerably nontrivial distribution 
on the actual phase space of constrained variables, which is a "distorted" manifold (or a 
leaf) immersed in the total space. Interesting structures created in a fluid or a plasma may 
be delineated by unearthing leaves of the phase space and analyzing their distortion with 
respect to the energy norm. When one can "integrate" the kernel of the Poisson operator 
to construct Casimir elements, the Casimir leafs foliate the Poisson manifold and, then, the 
effective energy is the energy-Casimir functional. 
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In addition, here we have proposed a model for a physical process that removes the con- 
straints of Casimir elements and enables the system to seek lower-energy states on different 
Casimir leaves. By invoking an extended phase-space, we canonized the Poisson operator 
and introduced a coupling of the original ideal system with an external energy source — the 
exchange of energy between the original system and the connected external system thus 
described a "dissipation" process. This formulation is based on the method of "minimum 
canonization" that interprets Casimir elements as "adiabatic invariants," and "unfreezes" 
the Casimir elements to be dynamic, by perturbing the Hamiltonian with respect to the new 
angle variable added to the phase space; such perturbations increase the number of degrees 
of freedom and are therefore a kind of singular perturbation. 

The theory was applied to the tearing-mode instability, where a tearing mode was re- 
garded as an equilibrium point on a helical-flux Casimir leaf. As long as the helical-flux is 
constrained, the tearing mode cannot grow. However, it was shown that a singular pertur- 
bation that allows the system to change the helical flux can cause a tearing mode to grow if 
it has an excess energy with respect to a fiducial energy of the Beltrami equilibrium at the 
bifurcation point. 
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